Astronomy & Astrophysics manuscript no. paper 


February 2, 2008 


(DOI: will be inserted by hand later) 





Catastrophic alpha quenching alleviated by helicity flux and shear 

Axel Brandenburg^ and Christer Sandin^ 



^ NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

^ Astrophysical Institute Potsdam, An der Sternwarte 16, D-14482 Potsdam, Germany 



o 
o 

(N 



0^ 
(N 

> 

(N 
O 

o 
6 



Received 16 January 2004 / accepted 21 



July 2004 



Abstract. A new simulation set-up is proposed for studying mean field dynamo action. The model combines the 
computational advantages of local cartesian geometry with the ability to include a shear profile that resembles the 
sun's differential rotation at low latitudes. It is shown that in a two-dimensional mean field model this geometry 
produces cyclic solutions with dynamo waves traveling away from the equator - as expected for a positive alpha 
effect in the northern hemisphere. In three dimensions with turbulence driven by a helical forcing function, an alpha 
effect is self-consistently generated in the presence of a finite imposed toroidal magnetic field. The results suggest 
that, due to a finite fiux of current helicity out of the domain, alpha quenching appears to be non-catastrophic 
- at least for intermediate values of the magnetic Reynolds number. For larger values of the magnetic Reynolds 
number, however, there is evidence for a reversal of the trend and that a may decrease with increasing magnetic 
Reynolds number. Control experiments with closed boundaries confirm that in the absence of a current helicity 
fiux, but with shear as before, alpha quenching is always catastrophic and alpha decreases inversely proportional 
to the magnetic Reynolds number. For solar parameters, our results suggest a current helicity fiux of about 
0.001 G^/s. This corresponds to a magnetic helicity flux, integrated over the northern hemisphere and over the 11 
year solar cycle, of about 10*^ Mx^. 
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1. Introduction 

The large scale magnetic field of stars and galaxies is often 
interpreted in terms of mean field dynamo theory, which 
takes into account that the turbulence is at least partially 
helical. The helicity, in turn, leads to the so-called a effect, 
i.e. an averaged field-aligned current that induces new field 
loops perpendicular to the original field (MofFatt 1978, 
Krause & Radler 1980). In the presence of shear, poloidal 
loops get sheared out, thereby reinforcing the toroidal field 
from which more poloidal loops can be created. This is the 
basic afl dynamo mechanism which is at least in principle 
able to explain the cyclic variation of the solar magnetic 
field (Parker 1979). 

The investigation of mean field dynamos has been 
the subject of numerous papers since the 1970s (see also 
Zeldovich et al. 1983). With the advent of high resolution 
turbulence simulations it has become exceedingly clear 
that there is a serious problem in the nonlinear case at 
large magnetic Reynolds numbers. Similar problems are 
completely unknown in the context of nonmagnetic, purely 
hydrodynamic turbulence or in nonhelical hydromagnetic 
turbulence where the associated dissipative fluxes always 
remain finite. This is not the case with the magnetic helic- 
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ity fiux which goes to zero in the large magnetic Reynolds 
number limit (Berger 1984). The magnetic helicity is thus 
almost perfectly conserved in practically all astrophysi- 
cally interesting cases. 

An important consequence of magnetic helicity conser- 
vation is the fact that the a effect cannot produce any net 
magnetic helicity. This means that if the a effect produces 
large scale magnetic fields, it must at the same time also 
give rise to a certain amount of small scale fields with op- 
posite sign of magnetic helicity (Seehafer 1996, Ji 1999). 
Hence the strength of the large scale field that can be gen- 
erated on dynamical time scales is limited, as the associ- 
ated small scale field cannot grow significantly above the 
equipartition field strength (Brandenburg 2001, hereafter 
BOl, Field & Blackman 2002, Blackman & Brandenburg 
2002, Subramanian 2002). Recent work has shown that 
this leads to a rather restrictive nonlinearity of the a ef- 
fect. It has been recognized for some time (Blackman & 
Field 2000a, b, Kleeorin et al. 2000) that the conservation 
of magnetic helicity may be particularly restrictive in the 
presence of closed or periodic boundary conditions used 
in many investigations. While it is clear that magnetic 
helicity fluxes through boundaries can help in principle 
(Brandenburg et al. 2002), one must still ensure that it 
also is of the right properties. If the magnetic helicity fiux 
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Fig. 1. A sketch of the solar angular velocity at low lat- 
itudes with spoke-like contours in the bulk of the convec- 
tion zone merging gradually into uniform rotation in the 
radiative interior. The low latitude region, modeled in this 
paper, is indicated by thick lines. 



carries away most of the desired large scale field, noth- 
ing will be gained and the dynamo will be even less effi- 
cient. This is indeed what early simulations have shown 
(Brandenburg & Dobler 2001). 

The purpose of the present paper is to show that the 
situation changes considerably when helicity flux is medi- 
ated by shear, the only mechanism known that can sep- 
arate and hence also transport magnetic helicity in space 
(the a effect, by comparison, separates and transports 
magnetic helicity in wavenumber space). 

We adopt the simplest possible model that is able to 
capture the effects of helicity transport from one hemi- 
sphere to the other. To motivate our model we first look 
at an idealized representation of the solar angular velocity 
which is spoke-like in the bulk of the convection zone and 
nearly rigid in the radiative interior; see Fig. ^ 

We model the region below 30° latitude by adopting a 
cartesian geometry where the x direction corresponds to 
radius, the y direction to longitude, and the z direction 
to latitude. We ignore the fact that in the sun the radial 
transition to uniform rotation is much steeper and model 
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Fig. 2. Differential rotation in our cartesian model, with 
the equator being at the bottom, the surface to the right, 
the bottom of the convection zone to the left and mid- 
latitudes at the top. The boundary conditions for the three 
components of the magnetic field and the vector potential 
(indicated near the boundaries of the box) are discussed 
in Sect. 121 

the mean toroidal velocity simply in terms of trigonomet- 
ric functions using 

U = Uocoskixcoskiz, (1) 

where ki is the lowest wavenumber in the (x, z) plane with 
— 7r/2 < kix < and < kiz < Tr/2. In the following we 
adopt units where ki = 1. The equator is assumed to be 
at z = and the outer surface at x — 0. The bottom 
of the convection zone is at a:; = — 7r/2 and the latitude 
where the surface angular velocity equals the value in the 
radiative interior is at z = 7r/2; see Fig.|21 

In order to clarify some basic properties of this rather 
unexplored geometry we begin by studying afl dynamos 
in this geometry in Sect.|21 after introducing our numer- 
ical approach in Sect. |21 The problem is at first reduced 
to two dimensions. Here we assume that there is an a 
effect so that a mean magnetic field can readily be gener- 
ated by this term. Next we drive helical turbulence by a 
corresponding forcing term in three dimensions and mea- 
sure the resulting a effect by imposing a toroidal magnetic 
field in Sect.0] In particular, we study the dependence of 
a quenching on the magnetic Reynolds number and show 
that the a effect is only catastrophically quenched when 
the boundaries are closed or when the large scale shear 
has been turned off. 

2. Our numerical approach 

The evolution of the magnetic field, B, is governed by the 
induction equation, dB/dt = — V x where E is the 
electric field. The induction equation is solved in terms of 
the magnetic vector potential. A, where B ~ "V x A. The 
evolution of A is therefore given by dA/dt = —E — V0, 
where (f> is the scalar potential. By using the gauge in 
which (j) ~ we simply have 

dA/dt = -E, (2) 
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which is the equation that will be considered throughout. 
The (negative) electric field is given by 

E = U X B 77^0 J, (3) 

where U is the velocity, J = 'VxB/ fiQ the current density, 
and fiQ the vacuum permeability. 

At the bottom of the convection zone and at mid- 
latitudes we assume the same boundary conditions for 
the magnetic field as for the velocity field, i.e. the field 
is tangent to the boundaries, which corresponds to per- 
fect conductor boundary conditions. This means that the 
tangential electric field vanishes, and, because of the gauge 
= 0, we have 

Ax = Ay = on z = 7r/2, and (4) 

Ay ^ A^ ^0 onx = -tt/2. (5) 

On the equator and at the outer surface we assume that 
the magnetic field is normal to the boundaries, i.e. 

Bx ^ By = on z — 0, and (6) 

By = B,^=0 on X = 0. (7) 

At the equator this boundary condition is consistent with 
dipolar parity of the field. The full set of boundary condi- 
tions for all three components of both A and B are given 
on the four 'meridional' boundaries of the box in Fig. [21 
In the y direction we adopt periodic boundary conditions 
over the interval < y < 2tt. 

Occasionally we refer to the boundary conditions (jSJ 
and 10 as open, because they permit flux of magnetic and 
current helicities through the z = and x ~ boundaries. 
For comparison we also perform calculations with closed 
boundaries where A^ = Ay = on z = and Ay = Az = 
on X = 0. With these conditions the magnetic and current 
helicity fluxes vanish through these boundaries. At the 
bottom of the convection zone and at mid-latitudes the 
fluxes of magnetic and current helicities are always van- 
ishing. This is probably a reasonable assumption, because 
at these locations there is no shear to mediate large scale 
helicity transport and the small scale helicity transport 
was already previously found to fluctuate around zero if 
there is no shear (Brandenburg & Dobler 2001). 

For both the mean field calculations and the turbu- 
lence simulations we step the equations forward in time 
by using the Pencil Codc"'^ . For the mean field calculations 
a typical resolution of 32^ meshpoints proved to be suffi- 
cient, while for the turbulence simulations the required 
resolution depends on the magnetic Reynolds number, 
i?in. Here, i?ni is based on the magnitude of the turbu- 
lent velocity and not on the larger shear flow velocity. 
(The precise definition is given below in Sect. |121) For 
i?in ~ 100, a resolution of 512'^ meshpoints is required, 
while for R^^ « 15, a resolution of 128"^ meshpoints proved 
to be sufficient. We note, however, that the aspect ratio 

^ [http : //www.nordita.dk/sof tware/pencil-codel 
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of the box is 1:4:1 and, although there is shear smear- 
ing out structures in the y direction, a uniform mesh as- 
pect ratio seems often to be preferred. For example, a run 
with 128 X 512 x 128 meshpoints allowed us to use higher 
Reynolds numbers than 128"^ meshpoints. 

3. Mean field calculations 

An important aspect of our studies is to show that the 
shear flow depicted in Fig.|21is a reasonable approximation 
to the differential rotation present in the sun (which is 
more like that depicted in Fig.^. 

In the context of mean field theory it is known that 
in spherical shells both dipolar and quadrupolar solutions 
are approximately equally easily excited (Roberts 1972) 
and that both solutions can be oscillatory with field mi- 
gration away from the midplane (when a > 0). We want 
to know whether in the present geometry the magnetic 
field evolution is similar to that in spherical shells. 

The mean field induction equation is given by dB /dt = 
— V X which again is solved in terms of the mean mag- 
netic vector potential, i.e. dA/dt = —E, where 

~E = TJxB + £- rjnoJ. (8) 

Here, overbars denote azimuthal averages, and velocity 
and magnetic field arc split into mean and fluctuating 
components via B = B + b and U = U + u. The electro- 
motive force from the fluctuating components of velocity 
and magnetic flcld, 5 = it x 6, is in its simplest form (e.g. 
MoflFatt 1978, Krause & Radler 1980) 

£ = aB - ritfioJ, (9) 

where a (related to the mean helicity) and rjt (turbulent 
diffusivity) could still be functions of x, z, and t, as well as 
a function of B itself, but for simplicity we assume them 
to be constant here. 

The solutions are characterized by two non- 
dimensional parameters, 

Ca = a/{T]Tki), and Cs = Uo/{riTki), (10) 

where r]^ = rj + rit is the sum of microscopic and turbulent 
magnetic diffusivities. 

In Fig. Owe plot the stability diagram in the {Ca, Cs) 
plane. For Ca < Ca,crit the solutions are decaying and 
for Ca > Ca^crit they are growing exponentially and are 
oscillatory (Hopf bifurcation) , except for a narrow interval 
around Cs = 0. Such a behavior is quite typical of aQ 
dynamos (see, e.g., Roberts & Stix 1972). For Cs = 1000 
we have also considered the quadrupolar solution and flnd 
that it is slightly easier to excite (see Fig. OJ. As stated 
earlier, the approximately equal excitation conditions for 
dipolar and quadrupolar solutions, seen in Fig.|31 is typical 
of aH. dynamos in spherical shells. Indeed, the fact that 
quadrupolar solutions can be preferred has been found in 
other solar dynamo models (Dikpati & Oilman 2001). 

In Fig. 01 we show contours of By for the marginally 
excited case with Cs = 1000 for different times covering 
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Fig. 3. Critical value of Cq for dynamo action as a func- 
tion of shear parameter, Cs- Note the typical decrease of 
the critical value of Cq with increasing Cs- 



a little more than half a cycle. One clearly sees magnetic 
field migration away from the equator. We note that in the 
sun the field migration is toward the equator. The reason 
for this is not entirely clear, but it could be caused by a 
negative a effect (but the reason for this is not clear ei- 
ther) or by meridional circulation (Choudhuri et al. 1995, 
Durney 1995). However, before these questions can seri- 
ously be addressed, it is mandatory to have a reliable mean 
field theory. This has so far been hampered by not being 
able to model the nonlinear feedback correctly. 

Based on the similarity of the stability diagram and the 
migration pattern in calculations using cartesian geometry 
and spherical shells (e.g. Roberts & Stix 1972), we may 
conclude that our local model provides a reasonable ap- 
proximation to the more realistic case of a spherical shell. 
Since azimuthal averages are used, the mean field equa- 
tions would be axisymmetric and could easily be solved 
using meshpoint methods. The case of three-dimensional 
turbulence in spherical shells is considerably more difh- 
cult because the coordinate singularity at the poles leads 
to serious timestep restrictions in the azimuthal direction. 



4. Turbulence simulations 

In this section we consider three-dimensional turbulence in 
the same cartesian geometry that we used in the previous 
section. We now consider the full (non-averaged) velocity 
and magnetic fields, U and B, respectively. 



4.1. Basic equations 

The fluctuating velocity together with the shear flow must 
be obtained by simultaneously solving the induction equa- 
tion [Eq. (O] together with the momentum equation. 
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Fig. 4. Contours of By for times separated by At — 300 
in the temporal sequence tl-t6. Note the field migration 
from the lower left corner to the upper right; in this figure 
C's = 1000 [see Eq. 



which we write here for an isothermal gas of constant 
sound speed Cg, 



BU 
~Dt 



-c^Vlnp 



J X B 

P 



f 



(11) 



where / is the forcing function driving both the turbulence 
(around a narrow band of wave numbers around fcf = 5) 
and the shear flow (cf. Brandenburg et al. 2001). Moreover, 



N^U + i VV • (7 + 2S • Vlnp) 



(12) 



is the viscous force where — const is the kinematic vis- 
cosity, Sij — \ [Uij + Uj^i) — ^SijUk^k the traceless rate of 
strain tensor, and p the density which obeys the continuity 
equation which we solve in the form 



Dlnp 
Bt 



-V-C/, 



(13) 



where D/Di = d / dt + U ■ 'V is the advective derivative. 
We adopt a random forcing function with finite helicity; 
see BOl for details. 
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4.2. Results for a quenching 

We have carried out a range of simulations for different 
values of the magnetic Reynolds number, 

Rm = u,„,J{r]kf), (14) 

for both open and closed boundary conditions. (Here, Urms 
does not include the mean shear flow.) In order to measure 
a, a uniform magnetic field, Bq = const, is imposed, and 
the magnetic field is now written as B = Bo+'VxA. In all 
cases presented below we have used Bq w O.OTyZ/IopUrms- 
In Fig. we show a graphical presentation of a typical 
snapshot of a solution. 

We have determined a by measuring the turbulent 
electromotive force, i.e. a = (S) ■ Bq/Bq. Similar inves- 
tigations have been done before both for forced turbu- 
lence (e.g., Cattaneo & Hughes 1996, see also BOl) and 
for convective turbulence (e.g., Brandenburg et al. 1990, 
Ossendrijver et al. 2001). The contribution from rj^J has 
been ignored in this approach; this is justified because for 
strong shear the poloidal field, giving rise to a toroidal 
mean current, is weak. 

It is well known that the a effect is an extremely 
noisy quantity - especially in the case of large magnetic 
Reynolds numbers (Cattaneo & Hughes 1996). The strong 
fiuctuations are also clear from Fig.|Bl where we plot a{t) 
both for open and closed boundaries. Note that a/urms 
fiuctuates in time in the range ±0.2 about a much smaller 
average value of, e.g., —0.03 in the top panel. 

As expected, a is negative when the helicity of the forc- 
ing is positive, and a changes sign when the helicity of the 
forcing changes sign. For ^ 14 the magnitudes of a 
begin to be different in the two cases: |a| is larger when 
the helicity of the forcing is negative. In the sun, this cor- 
responds to the sign of helicity in the northern hemisphere 
in the upper parts of the convection zone. This is here the 
relevant case, because the differential rotation pattern of 
our model also corresponds to the northern hemisphere. 

There is a striking difference between the cases with 
open and closed boundaries which becomes particularly 
clear when comparing the averaged values of a for differ- 
ent magnetic Reynolds numbers; see Fig. |7| With closed 
boundaries a tends to zero like R^^, while with open 
boundaries a shows no such immediate decline; only for 
larger values of i?m there is possibly an asymptotic a cc 
dependence. There is also a clear difference between 
the cases with and without shear. In the absence of shear 
(dotted line in Fig.d) a declines with increasing R^^ , even 
though for small values of R^ it is larger than with shear. 
This suggests that the presence of shear combined with 
open boundaries might be a crucial prerequisite of dy- 
namos that saturate on a dynamical time scale. 

The difference between open and closed boundaries 
will now be discussed in terms of a current helicity flux 
through the two open open boundaries of the domain. 
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Fig. 5. Vectors of U and B in an xy plane through 
z = 7r/4, superimposed on a grey/color representation of 
their normal components, for a run with R-^ = 14, open 
boundaries, shear, and negative helicity. Note that the ve- 
locity field is dominated by the toroidal shear flow. 

4.3. Current helicity flux 

It is suggestive to interpret the above results in terms 
of the dynamical a quenching model, where a is propor- 
tional to the difference between kinetic and current helic- 
ities (Pouquet et al. 1976), i.e. 

a = -\t [(jj ■ u - p^^j -b) = aK + au- (15) 

Here, a; = V x it is the small scale vorticity, j = V x b//io 
is the small scale current density, ax is the kinematic con- 
tribution to the a effect, and au the magnetic contribu- 
tion primarily responsible for the quenching of the a effect. 

In order to obtain an expression for we begin 
with the evolution equation for j ■ b. In the absence of 
boundary conditions it has been advantageous to start 
with the evolution equation for magnetic helicity (because 
it is conserved) instead of current helicity (which is not 
conserved); see Kleeorin & Ruzmaikin (1982), Klceorin 
et al. (1995), Blackman & Brandenburg (2002). In the 
case of open boundary conditions, this approach becomes 
cumbersome, because one has to consider the gauge- 
invariant relative magnetic helicity (Berger & Field 1984). 
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Fig. 6. Time traces of a{t) for runs with open boundaries 
(for both positive and negative kinetic hehcity) and closed 
boundaries (for positive kinetic helicity). i?ni — 30 in aU 
cases. 

Furthermore, the concept of a density of magnetic helic- 
ity is not meaningful, because it would depend on the 
gauge. In order to avoid these problems it is advantageous 
consider the current helicity equation (Subramanian & 
Brandenburg 2004, Brandenburg & Subramanian 2004). 
Apart from this technicality, the following derivation is 
similar to that of Blackman & Brandenburg (2002) for 
the case without current helicity flux. 

Using the evolution equation, db/dt = —V x e, for the 
fluctuating magnetic field, where e = E — E is the small 
scale electric field and E — j]J — £ the mean electric field, 
we can derive the equation 



where 

^SS 

c 



-2eTc-V 



C 1 



2e X ji + (V X e) X b 



(16) 



(17) 



is the current helicity flux from the small scale field, and 
c = V X J the curl of the small scale current density. 
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Fig. 7. Dependence of |(Q!)|/urms on R-a^ for open and 
closed boundaries. The case with open boundaries and 
negative helicity is shown as a dashed line. Note that for 
i?m ~ 30 the a effect is about 30 times smaller when the 
boundaries are closed. The dotted line gives the result with 
open boundaries but no shear. The vertical lines indicate 
the range obtained by calculating a using only the first 
and second half of the time interval. 



6, where fcf is 



2 = J J. In the isotropic case, e • c ; 
the typical wave number of the fluctuations, here assumed 
to be the forcing wavenumber. Ignoring the effect of the 
mean flow on £ [as is usually done; but see Krause & 
Radler (1980) or Rogachevskii & Kleeorin 2003], we obtain 

e~6 w -{u X Bn) ■ b + 7i]~b = £ B + rjj~b, (18) 



where we have used u x b = £ and Bq = B. Using 
standard expressions for the turbulent magnetic diffusiv- 
ity, rjt = irit^j^g, and the equipartition field strength, 



cq — \/ fJ-opUrms, wc eliminate r via 

Vt/B!^. (19) 
This leads to an explicitly time dependent formula for a, 
'I B 



da 
'di 



B^ 



(20) 



This equation is similar to that of Kleeorin et al. (2000, 
2002, 2003) who considered the flux of magnetic helicity 
instead of current helicity. 

Making the adiabatic approximation, i.e. putting the 
rhs of Eq. H20() to zero, one arrives at the algebraic steady 
state quenching formula {da/dt = 0) 



a 



aK + Rm {ritP-aJ ■ B 



1 + /B2 



(21) 



Furthermore, if the mean field is defined as an average 
over the whole box, then B = Bq = const, so J = and 
Eq. reduces to 

" ■ ^ " ' '^^'i (22) 



.v-y^^/Bi 



1 + RmBt/Bl 



oq 
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This expression applies to the present case, because we 
consider only the statistically steady state and we also 
define the mean field as a volume average. 

For closed boundaries, V • = 0, and so Eq. 
clearly reduces to a catastrophic quenching formula, i.e. a 
vanishes in the limit of large magnetic Reynolds numbers 
as 

ttK 



(closed) 



1 + 



(for i?m ^ oo). (23) 



eq 



The dependence suggested by Eq. is confirmed 
by the simulations (compare with the dash-dotted line in 
Fig. IT)). On the other hand, for open boundaries the limit 
i?,n — > oo gives 

(for ^ (»), (24) 



^SS 



-(V-JP; 



c )/(2fct'B^) 



which shows that losses of negative helicity, as observed 
in the northern hemisphere of the sun, would enhance a 
positive a effect (Kleeorin et al. 2000). In the simulations, 
the current helicity flux is found to be independent of 
the magnetic Reynolds number. This explains why the a 
effect no longer shows the catastrophic dependence 
(see Fig. El). 

4.4. Estimates for the Vishniac-Cho flux 

Theoretical estimates for magnetic helicity fluxes have 
been proposed by Kleeorin et al. (2000) and Vishniac & 
Cho (2001). The two fluxes are rather different. The ex- 
pression of Vishniac & Cho (2001) has been confirmed in- 
dependently and can be written in the form (Subramanian 
& Brandenburg 2004; see also the review by Brandenburg 
& Subramanian 2004) 



(25) 



In Fig. IHl we plot the profiles of this flux, averaged in the 
y direction, on the two open boundaries. We also show 
the time evolution of the averaged fluxes on the two open 
boundaries. It turns out that the magnitude of the two 
fluxes is large compared with = Wrmsfcf^oj tint the 
fluxes also fluctuate strongly in time, so it is important to 
average over long times. Furthermore, there is a clear ten- 
dency for the difference between incoming flux at the equa- 
tor (^^y*^, dashed line) and outgoing fluxes at outer sur- 
face (^y*^, solid line) to cancel partially, giving a smaller 
net flux. Nevertheless, since |^^^| > I^Y^I, the net out- 
going flux is negative, as expected for the northern hemi- 
sphere. 

4.5. Large scale current helicity flux 

In earlier work (Brandenburg & Dobler 2001) it was re- 
ported that the contribution to the magnetic helicity flux 
was outweighed by a much larger flux from the large scale 
field. In the present paper we work instead with the cur- 
rent helicity, and the current helicity flux from the large 
scale field is 




7r/8 7r/4 Stt/S 

k\z (or -k\x) 




1000 1500 

lirms ^ 1 i 



2000 



=vc. 



Fig. 8. The Vishniac-Cho flux on the outer surface {T^ ) 

and at the equator {T^ ), averaged over the y direction, 
for the run with a > and i?m — 14. In the upper panel 
the flux is also averaged over t and shown as a function 
of z and x, respectively. In the lower panel, the x and z 
components of the fluxes are averaged over the z and x 
directions, respectively, and shown as functions of t. The 
fat solid and dashed lines denote the running means for 
the two functions. 

For the vertical field condition, see Eqs ^ and O, the 
second term in Eq. H2t)|) vanishes. Assuming isotropy, the 
contribution from the first term involves £ — aB — rjtJ, 
but this does not contribute either, because aB x J does 
not have a normal component on the boundaries, and 
rjtJ X J = 0. In the first term only the mean fiow term 
contributes, so we get 



-2{U XB)XJ= -2{jyUy)B. 



(27) 



2E X J + {V X E) X B. 



(26) 



Inspection of the data suggests that this is indeed a good 
approximation and that therefore even in the simulation 
the normal component of 5 x J is nearly vanishing on the 
boundary; see Fig. |^ Both on the outer surface and on 
the equator is such that it corresponds to a loss of 
negative current helicity. 

In Fig. O we also show the small scale current helic- 
ity fluxes on the two boundaries (fat lines). There is a 
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Fig. 10. Sketch illustrating the directions of large scale 
(LS) and small scale (SS) negative current helicity fluxes 
and their approximate magnitudes (in units of ^o)- Note 
that at the outer surface negative current helicity is 
ejected both via small and large scale fields, while at the 
equator the contributions from small and large scale fields 
have opposite sign. The small scale losses at surface and 
equator partially cancel, giving a net loss of negative cur- 
rent helicity of only about 4 JFq. 



"^rms '''1 f 

Fig. 9. Normal components of the current helicity flux on 
the outer surface (Tx) and at the equator (J-'z), averaged 
over the corresponding surfaces, for the same run as in 
Fig-El The fat lines denote the fluxes from the small scale 
field, ^(^, while the thin lines denote the fluxes form the 
large scale field, '^c^- The dotted lines near the two 
curves show the result of the approximation H27|) . 



tendency for the difference between incoming flux at the 
equator (dotted line) and outgoing fluxes at outer surface 
(solid line) to cancel, but the net outgoing flux is again 
negative. The flux for the total field is approximately four 
times larger than what is accounted for by the Vishniac- 
Cho flux. This might indicate that there is either another 
contribution to the current helicity flux, or that the r in 
the Vishniac-Cho flux is underestimated. 

4.6. Application to the sun 

The purpose of this section is to put some real numbers 
into the expression for the current helicity flux. Our sim- 
ulations have shown that a reasonable estimate for the 
current helicity flux at the outer surface is 

«30u™sfcfB2. (28) 

Applying this to the sun using Umis ~ 50m/s for the rms 
velocity in the deeper parts of the convection zone, fcf ~ 
10~^ cm~^ based on the inverse mixing length, and Bg « 
3 G for the mean field at the solar surface, we have w 
10~^ G^/s. The current helicity flux integrated over the 
northern hemisphere of the sun is then 4x 10^^ G^ cm^ s^^. 
Integrated over the 11 yr solar cycle we have 10^* G^ cm^. 

For the sun only magnetic helicity fluxes have been de- 
termined. As a rough estimate we may use « k^^^^ 
for the magnetic helicity flux. Using the same estimate for 
k{ as above we obtain about 10'*® Mx^ over the 11 yr solar 
cycle. This is indeed comparable to the magnetic helic- 
ity fluxes estimated by Berger & Ruzmaikin (2000) and 
DeVore (2000). 



We emphasize that, in the present context, "small 
scale" means 27r/fcf « 60 Mm, i.e. about one pressure scale 
height at the bottom of the convection zone. By contrast, 
"large scale" refers to length scales of the order of several 
hundred Mm, which is typically beyond the scale captured 
in the usual vector magnetograms. 

5. Conclusions 

The present simulations give a clear indication that the 
proposed set-up with open boundaries and shear can al- 
leviate the catastrophic quenching problem in mean field 
dynamo theory. Already in the absence of shear the resul- 
ting a effect can be larger when open boundary conditions 
are used, but there a still decreases with increasing mag- 
netic Reynolds number. In the presence of shear, allevi- 
ated catastrophic quenching is associated with a net loss 
of small scale current helicity. At the equator, negative 
small scale current helicity flows into the northern hemi- 
sphere, but there is an even larger negative small scale 
current helicity flows ejected at the outer surface. This 
results in a net loss of negative current helicity from the 
northern hemisphere, and corresponds to an integrated 
magnetic helicity flux of about 10^® Mx^ over the 11 yr 
solar cycle. 

Although our results are certainly encouraging, they 
must still be considered preliminary. First of all, we have 
still only considered a relatively limited range of magnetic 
Reynolds numbers; larger values are necessary before one 
can tell whether or not an asymptotic _Rm dependence still 
develops for larger values of i?m- Secondly, turbulence sim- 
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ulations of a dynamo are necessary to show that the elec- 
tromotive force that here is interpreted in terms of an a 
effect is indeed capable of generating large scale magnetic 
field of the type shown in Sect. 13 Of particular importance 
is the question whether the dynamo is oscillatory (as ex- 
pected from our mean field calculations) and whether the 
cycle frequency is independent of i?m- Especially at large 
values of Rm the direct approach tends to be advantageous 
compared to calculating a in the presence of an imposed 
field, because the resulting a is always much more noisy 
than the actual mean field obtained in a simulation (BOl). 

Another aspect to keep in mind is the fact that the 
correct boundary conditions for the solar dynamo are cer- 
tainly more complicated than the vertical field condition 
adopted here. There are good reasons to believe that the 
sun loses significant amounts of magnetic and current he- 
licity via coronal mass ejections (DeVore 2001, Demoulin 
et al. 2002a, b, Gibson et al. 2002; see also Blackman & 
Brandenburg 2003). The losses via via coronal mass ejec- 
tions are not easy to model within the present approach. 
An obvious possibility is to include an outer layer that re- 
sembles some important aspects of the solar corona (low 
density and hence low plasma beta). It is as yet unclear to 
which extent a spherical geometry is important. The solar 
wind, for example, cannot be modeled in cartesian geom- 
etry, and in some sense coronal mass ejections are just a 
particularly bursty and localized manifestation of the so- 
lar wind. It may therefore be worthwhile to consider the 
effects of boundaries in global simulations. 
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